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CN I Abstract 

Form factor representation of the correlation function of the 2D Ising model on 
<^ I a cylinder is generalized to the case of arbitrary disposition of correlating spins. 

The magnetic susceptibility on a lattice, one of whose dimensions (N) is hnite, is 
■ calculated in both para- and ferromagnetic regions of parameters of the model. The 

singularity structure of the susceptibility in the complex temperature plane at finite 
values of N and the thermodynamic limit N ^ oo are discussed. 
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^ ■ 1 Introduction 

o . 

The Ising model has long been a subject of great interest. The computation of the free 
energy p^j and spontaneous magnetization [2], Toeplitz determinant representations [3], 
form factor expansions [1] and nonlinear differential equations O |6] for the correlation 
p • functions are among the most important advances of the modern mathematical physics. 
D . The partition function of the 2D Ising model in zero field was calculated exactly [7j 
^ ' not only in the thermodynamic limit but also for finite lattices with different boundary 
conditions. The simplicity of the corresponding expressions enables one to get an idea 
^ ' about the mechanism of appearance of critical singularities in thermodynamic quantities 
' from both mathematical and physical points of view. 

Analytical expressions for thermodynamic quantities, which contain the dependence 
on lattice size, have numerous applications. For example, in computer simulation of ther- 
modynamic systems or quantum field models one often needs such expressions to estimate 
the number of degrees of freedom for which a discrete numerical model is adequate to the 
initial continuous and infinite system. It is worth mentioning that modern experiments 
and technologies often deal with finite-size systems. Theoretical analysis of such systems 
experiences the lack of exactly solvable examples. 

In this paper we present exact expressions for the 2-point correlation function and the 
susceptibility of the 2D Ising model on a lattice with one finite {N = const) and the other 
infinite (M oo) dimension. These expressions are very similar to well-known form 
factor expansions [8], [9]. We investigate the singularity structure of the susceptibility for 
finite values of N and discuss the thermodynamic limit N ^ oo. 
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2 Correlation function ((t(0, 0)cr(a:, 0)) 

The Ising model on M x square lattice (Fig. 1) is defined by the hamiltonian H[o] 
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Figure 1: The numbering of the lattice sites and the variants of the disposition of corre- 
lating spins: a) along the cylinder axis, b) arbitrary disposition of spins on the lattice. 



r 

where the two-dimensional vector r = {x,y) labels the lattice sites: x = 1,2, ...,M, 
y = 1,2, . . . , N; the Ising spin cr(r) in each site takes on the values ±1; the parameter 
J > defines the energy of the coupling of adjacent spins. Shift operators Va;, act as 
follows: 

V^a(x, y) = a{x + l,y), Vya{x, y) = a{x, y + 1). 
Partition function at the temperature 

Z = ^e-^^M, (1) 

and 2-point correlation function 

{a{r,)a{r,)) = Z'' e-^'^^MriMr^) (2) 

are given by the sums over all possible spin configurations. It is convenient to introduce 
the following dimensionless parameters: 

K = (3J, t = tanhK, s = sinh2K. (3) 
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We will impose periodic boundary conditions on both axes. This gives two equations 
for the shift operators V^;, V^: 

(V.)^ = l, (V,)^ = l. 

For such boundary conditions the partition function ([T]) can be written as a sum of four 
terms [7]: 

Z = {2 cosh^ Kf''' . i {q^f'f^ + Q(^'^) + g(^'/) - Q(^'^)^ . (4) 

Each of them is given by the pfaffian of the operator D (a lattice analogue of the Dirac 
operator) 

Q = Pf5, (5) 

defined by 

/ 1 1 + 1 \ 

-1 1 1+tV 



D 



y 

-1 - tV-r -10 1 



(6) 



\ -1 -l-tV-y -1 / 

The upper indices (/, b) of the quantities Q in (jl]) correspond to different types (antiperi- 
odic or periodic) of boundary conditions for the operators V^^, Vy in 

iyfr = i^fY = 1, i^i'Y = (vf )^ = -1. (7) 

When, for example, M ^ N (i.e. the torus transforms into a cylinder), then in the right 
hand side of (jlj) only the "antiperiodic" term survives: 

Z = (2cosh2ir)^^^Q(/'/). (8) 

Since the operator D is translationally invariant, the pfaffian ([5]) can be easily computed. 
Using Fourier transformation, one finds the following factorized representation for the 
partition function (jS]): 

Z = 2^^^J]^^'^\s2 + l-s-cosq,-s- cosqyf'\ (9) 
q 

The superscript (/) in the products (or sums below) implies that the quasimomentum 
components and qy run in the Brillouin zone over half- integer values in the units 2tt/M 
and 27r/A^, respectively; integer values correspond to the superscript {h). For example, 

qy l=\ ^ ^ 9h l=\ ^ 

The product over one of the quasimomentum components in the right hand side of ^ 
can be calculated in an explicit form, so that for the partition function one has 

Z = (2s)^^/2-Q(/)g-A./,(,)/2(^ ^ e-*^^^'')), (10) 



where the function 7(g) is the positive root of the equation 



. ,2 7(g) -12/^, . 2 9 

smh ^ = smh — + sm -, 



(11) 



and the parameter /x is the following function of s: 

sinh ^ = -^(\/s — l/Vs 




(12) 



For q 7(g) remains positive in the whole range of the parameter s {0 < s < 00), but 
7(0) changes its sign after crossing the critical point s = 1. Since the product in flTUl) 
is taken over fermionic spectrum, which does not contain the value q = 0, this does not 
cause any problem there. However, we will see that the ambiguity in the definition of 
7(0) = ±/i leads to two different representations for the correlation function. 

The sum over spin configurations in the right hand side of the expression ([2]) for 
the correlation function can also be written in terms of pfaffians [lOj. Corresponding 
matrices, however, are not translationally invariant. This fact drastically complicates the 
calculations. Nevertheless, the computation of the correlation function can be reduced to 
the computation of the determinant of a matrix 



of considerably smaller size dim = dim(r), defined by the distance between correlating 
spins. Further work is needed to transform the representation (fT3|) into a representation 
with analytic dependence on the distance. 

Form factor representation for the correlation function of the Ising model is the most 
acceptable from physical point of view. First it was obtained in ^ for the infinite lattice in 
the ferromagnetic region {K > Kc, s > 1). Later it was extended [5] to the paramagnetic 
case {K < Kc, s < 1). We note that somewhat earlier a similar representation for the 
2-point Green function was deduced in [11] via S'-matrix approach [12] for a quantum 
field theory model with factorized S'-matrix (5*2 = —1), which is usually associated with 
the scaling limit of the Ising model. The discovery of the form factor representation for 
the correlation function has led to the whole trend [13] in the integrable quantum field 
theory. 

For the finite lattice the problem seems to be more difficult, but the result [H] is, in 
a sense, even simpler. If correlating spins are located along one of the lattice axes, the 
matrix in the right hand side of (fT3l) has Toeplitz form. For example, when correlating 
spins are located along the horizontal axis (Fig. la), then 



(o-(O)a(r)) = detA^'^''^\ 



(13) 



(a(ri)a(r2)) =detA(l^l), 



- ri 



(x,0) 



(14) 



and the elements of |x| x |x| matrix Ajl^j^, are given by [H] 

.{\x\) _ 1 ^UJ) e*P-('^-'=')[2t(l + t^) - (1 - t^){e^P- + t'^e-'P-)] 




(15) 



MN^ (1 + t2)2 _ 2t(l - t2)(cosp^ + cospy) 

J k —— • • • 7 1*^1 1* 
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As it was shown in [H] via Wiener- Hopf integral equations technique adjusted to 
finite-size lattice, the determinant f|T^ can be computed analytically and for the correla- 
tion function one obtains 

[N/2\ 

{a{v,)a{v,)) = (e-eT)e-'-'/^5^^72K^), for 7(0) = /z, (16) 

1=0 

(a(ri)a(r2)) = (^ eT)e-l^l/^ Yl 92i+i{x), for 7(0) = -/i, (17) 

1=0 

■E"{n'-^yM «o^i, (18) 



[1] 



FM = ni!fc*)(|, = (19) 
AAsmh((7i + 7j-)/2) 

where ji = j{qi), rji = r]{qi). The expressions ( fT6l) . ( fT7|l are finite sums. However, the 
upper limits of summation can be set infinite, since it follows from ( |T9l) that the form 
factor Fn[q] vanishes for n > N. Note an important detail — the summation over the 
phase volume in (ITSl) is taken over bosonic spectrum of quasimomenta, in contrast with 
the initial fermionic spectrum, which determines the matrix (fT5l) . The other quantities in 
(fT6ll -( fT9il are given by 

^ = |i_s-4|V', (20) 

TT 

sin((p + q)/2) 



^ f dpdqi{p)i{q) 

27r2 y sinh(A^7(p))sinh(iV7(g)) 



(21) 



sin((p- g)/2) 

A-^ = - I dplncoth{Nj{p)/2), (22) 



TT 







Viq) = - [ ^^^^^^f-^^\ncoth{N^{p)/2). (23) 
Tc J cosh7(gj — cosp 


"Cylindrical parameters" ^t, A~^, ri{q) explicitly depend on the number of sites on 
the base of the cylinder. Their asymptotic behaviour for A^|yu| ^ 1 is the following: 

InCr ^ -e-2^l^l, (24) 

TT 



A-1 ^ ^-^H |2sinh|/.| 



/ N 4e ^1^1 /sinh lul 

Thus outside the critical point cylindrical parameters A~^, In^^ and 77(g) exponentially 
decrease for large A^ and tend to zero for infinite lattice. Finite sums (1161) . (1171) transform 
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into series, summation over the phase volume in fllSI) is substituted by integration and, 
as a result, form factor representations on the cylinder transform into form factor repre- 
sentations on the infinite lattice [8], [9j. For any finite N both expansions — over even 
n ( fT6l) and over odd n ( fTTl) — are valid in both ferromagnetic (s > 1) and paramagnetic 
(s < 1) regions. However, for oo, the first series is well-defined and converges in the 
ferromagnetic region and the second one does so in the paramagnetic region. 

Recall that we started from the determinant f|T4|) of a |x| x |x| matrix. The number of 
terms in its formal definition rapidly increases when x grows. However, the form factor 
representations (fT6l) -(fT9l) are finite sums for any fixed A^, and the number of terms in 
these sums does not depend on This gives a unique opportunity to verify (|T6l) -( fT9l) 
by comparing these representations with the results of transfer matrix calculations for 
N-Tow Ising chains. For fixed N the dimension of the corresponding transfer matrix is 
equal to 2^ x 2^. One can find analytically all eigenvectors and eigenvalues if is not 
too large. We have successfully performed such check analytically for N = 2,3,4 and 
numerically — for = 5, 6. 

3 Correlation function ((j(0, 0)(j(x, y)) 

The rigorous derivation of the form factor representation on the cylinder was performed 
in [H] only for the spins located along the cylinder axis. We have not yet succeeded 
in generalization of the method for arbitrary disposition of correlating spins (Fig. lb). 
Meanwhile, the calculation of the momentum representation of correlation function 

G(p) = 5^e^P^(a(0)a(r)), (27) 

r 

or the susceptibility (which is related to G(p = 0)) requires an explicit dependence on 
both components of the vector r. Form factor representations f|T6l) - f|T9l) have a trans- 
parent physical content. This allows to make reasonable assumptions for corresponding 
generalizations. The above mentioned possibility of independent check allows to eliminate 
wrong hypotheses and to make correct choice. In principle, when ?/-component of the vec- 
tor r is not equal to zero, all quantities in ( fT6i) -(|T9 l) could change their form. However, 
corresponding expressions for free bosons and fermions on the lattice prompt one of the 
simplest generalizations — just the substitution 

g-|x|7(g) _^ Q-\x\^{q)-iyq 

Suprisingly enough, this turns out to be sufficient. If instead of gn{x) (fT8|) one uses the 
expression 

then the correlation functions (fT6|) and (fT7|) exactly coincide with the transfer matrix 
results for A^ = 2, 3, 4 in the whole range of the variables x, y, K. Numerical calculations 



/p-Fl7j-»y93-%- \ 
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confirm this for = 5, 6 as well. Tlie validity of (1281) is out of doubts and we hope that 
the known answer will simplify the problem of its rigorous derivation. 

As an illustration, consider the example of = 3. The expansions f|T6l) - f|T7j) are very 
similar to the representation of the correlation function in terms of the transfer matrix 
eigenvalues 

(a(O)a(r)) = ai(y)(Ai/Ao)l^'l + a2{y){\2/\,t\ + ■ ■ • , (29) 

where Aq is the largest eigenvalue and the coefficients aj{y) are given by some bilinear 
combinations of the components of eigenvectors. To reduce, for example, f|T7|) to fl29|) . we 
use the following expressions for the cylindrical parameters v{q)'- 



q Q 



(30) 



J]W sinh (^iM+lM^ 



n^^^ sinh 

9 



£4 _ 1 p ^ (32) 

n^^^n^'^ sinh (^(^) n^^^n^^^ sinh (^M+iM) ■ 

q p ^ ^ q p ^ ' 

One can derive these expressions from ( !2TI) - (!23l) by passing to contour integrals in the 
variable z = e*^ and computing the residues. 
For A^ = 3 we have from (l3nD-(l32D and (EOD 



A-i = i[7(7r) + 27(7r/3)-7(0)- 27(27r/3)], (33) 

sinh 2M±I(V3) sinh 7(Z0±I(M3) ^-^^2 7(2V3)+7(V3) 

ffr = (34) 

sinh I^toM) sinh sinh 7(7r/3) sinh 7(277/3) ' 

dnh 7(0)+7(g) „;„V,2 7(27r/3)+7(g) 

A-^-n{q) ^ smli ^ smn ^ 

sinh ^(^)+^('') sinh^ 7(^/3)+7(g) ' ^ ' 



Then one finds 



ln(Ao/Ai) = A-i+7(0), (36) 
ln(Ao/A2) = A-i+7(27r/3), (37) 
ln(Ao/A3) = A-i+7(0) + 27(27r/3), (38) 
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1 sinh sinh sinh^ 7(2V3)+7(./3) 

~ 3 sinh sinh sinh 7(71/8) sinh 7(27r/3) ' 

2 sinh ^(°)+^^(-/^) sinh ^^^^ _ , 

"^^^^ ^ 3 s.nh7(V3)s.nhli^ ----^^-^/^^' ^'^^ 
1 1 
= 64 sinh ^(°)+^^(-/^) sinh ^(-)+^(-/^) sinh ^W+f-/^) sinh ^W+f-/^) ^^^^ 



X 



sinh7(7r/3) sinh7(27r/3) sinh^ 7(V3)+7(2V3) " 



2 

Our 2^ X 2^ transfer matrix has 8 eigenvalues, and some of them are equal. Besides that, 
some eigenvectors have zero components. As a result, the expression for the correlation 
function (12^ contains only three (not seven) independent terms. If we take into account 
the definition ffTTl) . (|T2l) of the function 7(g) for particular values of quasimomentum 
q = 0, n/3, 27r/3, tt, we get an exact correspondence between these three terms and 

m-m. 



4 Momentum representation of the correlation 
function 

Since we have the expression fl28|) for 5'„(r), which explicitly depends on both components 
of r, we can make the Fourier transform. Let us write the momentum representation of 
([27]) in a form similar to ([TBD- ffT711 : 

G(p)=aTj]^n(p), (42) 
n 

5^.(p) = J]e-l-l/V(r)e^P^ (43) 

r 

where 

00 N 

E=EE- w 

r x=—oo y=l 

Performing the summation in fH3|) . we find 



sinhfA-i+f:7,VnM 



[q] S=i coshj^A-1 + ^ 7jJ - cosp^ ^ i=i 



The x-component of the quasimomentum has a continuous spectrum, G [— 7r,7r], but 
Py is discrete: 

= ^, 1=1,2 ...N. 
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Corresponding (5-function in the right hand side of fl45p has the meaning of the Kronecker 
symbol 



mod 



The function '(jrSs*) is periodic in p^., "Py with the period 27r. Inserting the "unity" 

A^^+n7(7r) 



1= y dujb^-^ ^^^j-^. 



A-i+n7(0) 

into the sum fH5l) (here 5 denotes Dirac 5-function) and interchanging the order of sum- 
mation and integration, we obtain 



A~l+n7(7r) 



-n/A 



Pn[UJ,Py) 



9n[P) 



rE'"(n 



sinhcj 

do; — pn{uj,py), 

COShCiJ — COSPa; 



(46) 



A-i+n7(0) 



F'nW(A-' + Y.^,-uj]6ipy-J2<lj)- (47) 
j=i / V j=i 



, sinh7,- 

[q] ^j = l 'J 

On the infinite lattice in the scaling limit the rotational symmetry is restored and fj46l) 
transform into the classical Lehmann representation in the quantum field theory. 



5 Magnetic susceptibility 

On M X N square lattice with equal horizontal and vertical coupling parameters the 
partition function Z depends on four variables 

^ ^ -f3H\cr]+hT,cr(r) 

Z = Z{K,h,N,M) = J2(^ ^ , (48) 

where dimensionless parameter h = (3T-C, Ti — magnetic field. The magnetization 971 and 
magnetic susceptibility x can be expressed in terms of the derivatives of the partition 
function with respect to h: 

aJl(A-,/.,iV,A/) = -ji^^ = (a>, (49) 

p-\(K,h,N,M) = ^ = J]({ff(0)ff(r)> - {af\ (50) 

The magnetization for /i = and finite M, N is equal to zero due to Z2-symmetry of the 
Ising model. This holds even when one of the dimensions is set infinite. In the last case, 
when, for example, M oo, N = const, 2D Ising model transforms into a ID chain with 
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N rows, for which the spontaneous symmetry breaking is impossible. The susceptibihty 
can be easily computed from p2l) - fH5l) 



X 
X 



[N/2] 

Xo+^ X21 

1=1 

^^tN coth(l/2A), 

[{JV-l)/2] 

X21+1 

1=0 

-n/A 



for 7(0) = /i, 



for 7(0) = -/i, 



(51) 
(52) 
(53) 



tE' ' n 



[9] 



sinh 7,- 

i=l " 



F^[q] coth 



i=l 



i=l 



In the paramagnetic region (s < 1) the expression (!53l) admits the limit N ^ 00 and 
tends to the susceptibility on the infinite lattice. However, in the ferromagnetic region 
{s > 1) one can consider the limit N —>■ 00 only for the quantity 



Xf = X-Xo 



00 

^X2h 
1=1 



(55) 



which reproduces well-known zero-field ferromagnetic susceptibility of the Ising model in 
thermodynamic limit. For large but finite the main contribution to the susceptibility 
is given by the term xo 



y sinh 



(56) 



which exponentially increases with the growth of the size of the cylinder base. It follows 
from (!56l) that the larger N — the smaller field 5h ~ e~^''^' is needed to order all spins 
on the lattice. 

Unfortunately, the exact solution for the partition function of the Ising model in 
external field is not known. However, the appearance of spontaneous magnetization can 
be deduced from the analysis of high- and low-temperature expansions. The rigorous 
definition of spontaneous magnetization is given by the following order of limits according 
to the Bogolyubov concept of quasiaverages: 



TloiK) = lim \ lim m(K, h, N, M)l . 



(57) 



However, if we conjecture the decreasing of correlations at large distances and the possibil- 
ity of interchanging of the corresponding limits, we can find exact solution for the squared 
spontaneous magnetization. It is equal to spin-spin correlation function fl20l) with infinite 
distance between correlating spins 



ml{K) = hm (a(O)a(r)) = (a(0))(a(oo)) = (a)' = ^. 

|r|^oo 



(5^ 
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Meanwhile, the sums over lattice of each summand in the right hand side of flSOl) do not 
converge in the thermodynamic limit. Therefore, the substitution of DJl^{K, 0, oo, oo) by 
the limiting value of the correlation function (which equals ^) under the (infinite) sum in 
the last step of the limits /i — > 0, M, N oo requires not only ( l58l) . but also the existence 
of the limit 

limj^^hm MN[m^{K, h, M,N)-^]} = f{K), (59) 



and, moreover, 



fiK) = 0. (60) 

The explicit dependence of the correlation function on the size A^, namely, the exponential 
tending of cylindrical parameters to their limiting values fl2^ - fl2^ . can be viewed as an 
argument in favor of the equalities flS^ . fl^ . 

The behaviour of the correlation function at large distances in the ferromagnetic region 
is mainly determined by the first term in the expansion (fT6|l . Note that it does not depend 
on y-projection of r 

Go(|r|)=aTe-l-l/^. (61) 

Therefore, the distance ~ A, for which spins are strongly correlated, rapidly increases (cf. 
with (125!) ) with the growth of A^. Physically it means that for "ferromagnetic" temper- 
atures the cylinder is divided into "domains" of size ~ A with nonzero magnetization, 
the magnetization of the whole infinite cylinder being equal to zero. It is clear that the 
squared spontaneous magnetization would be more naturally defined by the value of the 
correlation function at large distances |r| = R{N), which do not exceed the size of the 
domain 

N < R{N) < A. 

It follows from (!25|) that for sufficiently large these inequalities can be satisfied. In 
accordance with this, the sum over x with infinite limits in the definition of the thermo- 
dynamic limit of the susceptibility (1501) should be substituted by a sum with the limits 
that do not exceed the size of the domain. In this case the condition 

R N 
x=—R y=l 

can be treated as a formal substantiation of the definition (l55l) of the susceptibility in 
the ferromagnetic phase. We can now estimate the "thermodynamic cutoff parameter" 
RiN): 

R{N) < yXTiVe ~ e^l^l/2[7r/(2Arsinh |/i|)]'/^ 

We believe that these estimates slightly clarify the physical content of the formal thermo- 
dynamic limit procedure. 



6 Singularity structure 

The initial expression ([1]) for the partition function of the Ising model is a polynomial in 
s, and the solution is the factorized form of this polynomial. It provides an example 
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of the mechanism of Lee- Yang "zeros" [IS], which stipulates the appearance of critical 
singularities in the thermodynamic limit. The roots of the polynomial iQ are located on 
the unit circle |s| = 1 in the complex s plane. For any finite M and the zero s = 1 
on the real axis does not appear, since the fermionic spectrum does not contain the value 
of quasimomentum = Qy = 0. When one of the dimensions increases then zeros are 
concentrated on the circle \s\ = 1, forming a dense set. In the limit M — oo, = const 
they are transformed into a finite number (equal to N) of the root type branchpoints, 
located on the circle |s| = 1. To see this, one has to use the representation (fTOj) and the 
definition (fTTl) . (fT2l) of the function 7(g). These branchpoints, in turn, form a dense set 
with the growth of A^, but in the limit N ^ 00 they are transformed into four isolated 
logarithmic branchpoints s = ±1, ±i. As a result, the specific heat in the thermodynamic 
limit acquires a logarithmic divergence ~ In |1 — s|. It is worth noticing that the specific 
heat is expressed through the same function in both ferromagnetic and paramagnetic 
regions of s. 

One could think that a similar picture holds for susceptibility. Indeed, the initial ex- 
pression ([2]) for the correlation function for finite M and A^ is a ratio of polynomials in s. 
The formation of the singularities of the partition function, which stands in the denom- 
inator, we have just briefly described. Unfortunately, the polynomial in the numerator 
cannot be written in such simple factorized form. Nevertheless, our form factor represen- 
tation for M —>■ 00 and finite A^ shows that the correlation function has a finite number 
of root branchpoints on the circle |s| = 1. Their number is doubled in comparison with 
the case of partition function, since the expressions (fT6l) - (|T9l) . (1281) contain the functions 
7(9) f|TT|) - corresponding to both bosonic and fermionic values of quasimomentum. The 
susceptibility on the cylinder is given by the infinite sum of correlation functions and this 
can lead to the appearance of additional singularities. One can show, however, that these 
singularities do not appear on the first sheet of the appropriate Riemann surface. 

As an example, let us write down the susceptibility x dSI]) for A^ = 3, using the 
expressions (!39!)- (HT!) and representations (I33|) - (l35!) for cylindrical parameters 



7(0)+7(27r/3) -i^ 7(^)+7(27r/3) • , 2 7(27r/3)+7(7r/3) 

P-'X = 1^ . ' . ' ' coth All+IM + 



sinh ^^°)+^^"/') sinh ^^^^^fl^^ sinh7(7r/3) sinh7(27r/3) 

1 1 
+ 64 sinh sinh ^^^^^fl^) sinh ^W+f -/^) sinh -Z^) ^^^^ 

X ^ . i /3Wf2 /S^ COth( -^"^+^W+^^(^-/^) 

sinh7(7r/3) sinh7(27r/3) sinh' 7(^/3)+7(2^/3) V 2 



The singularities in s could appear due to zero denominator in (1621) . It is easily seen, 
however, that the corresponding factors 

. 7(9) +7(9') / / w . , 7(?) - 7(9') 
smh ^ = (cos q — cos q) j smh ^ 

for q ^ q' are not equal to zero. It can also be shown that on the first sheet of Riemann 
surface (which is determined by the condition of positivity of 7(13'), treated as functions 
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of s, for real s > 0) the arguments of cotangents in (1^2]) also have non-zero values: these 
factors appear as the result of summation over coordinate x. Therefore, the complete 
set of singularities of the susceptibility is exhausted by the branchpoints contained in the 
functions 



:,7(q) 



-(s + s-i) + sin^ ^ + \l^{s + s-i) - cos2 | 



(63) 



For each value of quasimomentum g 7^ 0, vr the function ( l63l) has four branchpoints. If we 



denote them by Sc 



Sr. e 



then 



cos (fc 



cos^ q/2 
- sin^ q/2 



(64) 



It is seen from ( !63|l . that for g = 0, vr there exist only two branchpoints Sc = ±i. One can 
now show that for any fixed N the total number of singularities is equal to AN — 2, and all 
singularities are located on the unit circle \s\ = 1. We represent the corresponding picture 
for = 3 in the Fig. 2. We do not discuss the limit N 00, when the singularities 





Figure 2: The location of the singularities of susceptibility x the complex plane s 
sinh 2/3 J for = 3. 



on the circle |s| = 1 form a dense set. This problem was seriously analyzed in [16] |17j. 
where the authors conjecture that the singularities form a natural boundary |s| = 1 for 
the susceptibility. 
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